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ABSTRACT 


This  work  is  intended  to  provide  automated  methodology  for  processing  seismic  and  infrasound  data  from 
scismo-acoustic  arrays  and  to  apply  the  methodology  to  regional  networks  for  validation  with  ground  truth 
information.  As  reported  last  year,  the  project  has  developed  automated  techniques  for  detecting,  associating  and 
locating  infrasound  signals.  Work  reported  this  year  focuses  on  refinement  of  both  the  detection  and  location 
component  of  the  program.  Ground  truth  information  developed  on  the  Korean  Peninsula  is  further  used  to  refine  the 
methodology. 

New  location  developments  arc  underway.  Wc  arc  incorporating  an  infrasonie  grid-search  location  method  into  a 
Bayesian  framework  that  provides  conditional  probabilities  on  the  location  given  the  observed  data  (intcr-array  time 
delays  and  backazimuths)  and  a  prior  distribution.  The  prior  provides  ancillary  information  on  the  event  location, 
which  can  be  as  simple  as  empirical  bounds  on  baekazimuth  deviations  and  group  velocities.  However,  wc  are 
investigating  extending  this  simple  approach  through  the  use  of  a  3D  ray-tracing  algorithm  to  impose  more  specific 
constraints  based  on  available  atmospheric  specifications. 

Finally  a  set  of  ground  truth  scismo-acoustic  events  have  been  developed  for  the  Korean  Peninsula  and  include  a 
large  number  of  events  from  a  limestone  quarry  in  the  Republic  of  Korea.  In-minc  seismic  and  acoustic  observations 
provide  the  basis  of  the  detailed  ground  truth.  These  events  occur  daily  with  multiple  events  on  active  days  and  span 
nearly  two  years,  providing  the  basis  of  assessing  spatial  and  temporal  variations  in  propagation  path,  including 
variations  in  baekazimuth  as  well  as  phase  and  group  velocity.  The  reported  analyses  of  these  data  provide  empirical 
bounding  constraints  for  the  location  algorithm  that  is  under  development. 
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OBJECTIVES 

There  are  two  primary  objectives  to  our  effort  on  the  second  year  of  this  contract.  Each  objective  builds  on  the 
automatic  algorithms  for  detection,  association,  and  location  of  infrasound  events  developed  in  the  first  year  of  this 
project. 

•  Development  of  an  infrasound  location  algorithm  with  a  firm  statistical  basis.  This  objective  is  a  natural 
extension  to  the  bounding-constraint  algorithm  developed  in  Year  1.  We  formulate  a  scheme  for  calculating 
likelihood  values  at  each  node  in  the  grid  search  method  based  on  the  agreement  between  observed  and  predicted 
values  for  backazimuth  and  arrival  time.  This  scheme  provides  the  opportunity  to  compute  robust  confidence 
bounds  on  infrasound  locations. 

•  Preliminary  analysis  of  infrasound  propagation  effects  using  ground-truth  data  from  Korea.  The 

groundt  ruth  dataset  being  acquired  for  this  project  provides  an  opportunity  to  test  and  validate  new  automated 
algorithms.  We  report  on  the  preliminary  analysis  of  this  dataset,  highlighting  the  effect  of  seasonal  variations  in 
atmospheric  winds  on  travel-times  in  the  Korean  peninsula.  Propagation  path  variations  such  as  those  documented 
in  this  study  can  be  incorporated  into  the  location  algorithm. 


RESEARCH  ACCOMPLISHED 


Infrasound  Location  Algorithm 

We  arc  currently  extending  the  grid-search  based  infrasound  location  methodology  of  Arrowsmith  et  al.  (2008a), 
discussed  in  last  years  proceedings  (Arrowsmith  et  al.,  2008b),  by  incorporating  it  into  a  Bayesian  framework.  This 
development  enables  us  to  compute  non-ellipsoidal  confidence  bounds  on  the  event  location  using  two  constraints: 
arrival  times  t=[w„]  and  backazimuths  6—  [#lv. 


Bayes’  theorem  is  given  by: 


P(t  x  V  vlt  f)\-  ’V)  p(f  \ 

*  0^  *  Vo>Wo»v;- 


(1) 


First,  we  formulate  the  value  of  the  likelihood  function,  P(t,(?j/0,.x0,jp0,v^  at  each  point  in  our  parameter  space. 
We  define  our  likelihood  function  as: 

^(t,^o>  Vo>v)=  n©M'o.  Vn.v)n°M^>*  Wo>v)  (2) 
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are  the  individual  likelihoods  for  the  Ith  array,  based  on  backazimuth  and  arrival  time  constraints  respectively  .The 
residual  terms,  y.  and  £*.,  are  given  by: 
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e,  =  t,  - 


.  ,  +(y,-y0) 

ln  T 


(6) 


where  the  fh  array  is  located  at  (x p^,),  (x0,y0J{i)  rePresenls  a  point  and  time  in  our  parameter  spaee,  and  V 
represents  the  group  velocity  that  is  assumed  to  be  constant  at  all  arrays  in  the  likelihood  term. 

The  uncertainties  in  baekazimuth,  CJ 0  and  <TV,  represent  the  standard  deviation  in  the  measurement  error  (a  function 
of  array  geometry  and  F-K  approach)  and  the  standard  deviation  in  the  error  produced  by  the  wind  respectively. 
Similarly,  the  uncertainty  in  arrival  time  estimation  is  denoted  by  <7^. 

We  impose  physical  constraints  on  the  model  parameters  in  the  Bayesian  prior  term,  P(jo9XOJyQ^V^  as  follows: 

|  positive  constant  0.22  <  v  <  0.34 

P  =  \  .  (7) 

0  otherwise 


This  constraint  assures  that  the  group  velocity  associated  with  our  solution  (x0,^0,/0)  lies  within  the  infrasonie 
range  (for  direct  waves  through  to  thermospheric  returns). 

The  final  solution  (or  posterior),  P(/0,JC0,J,0,  is  the  likelihood  of  the  event  location  and  origin  time  given 

our  observed  data  (arrival  times  and  baekazimuths).  This  likelihood  is  a  function  that  spans  some  pre-defined 
geographic  region  and  time  window.  We  apply  a  grid-seareh  method  in  order  to  find  the  location  and  origin  time 
associated  with  the  maximum  likelihood,  with  eonfidenee  regions  defined  by  the  likelihood  distribution 

Example 

We  apply  our  Infrasonie  Bayesian  Locator  (IBL)  to  a  ground-truth  roeket  motor  explosion  recorded  at  three  arrays  in 
Utah  (Stump  et  al.,  2007).  This  event  is  diseussed  in  Arrowsmith  et  al.  (2008a)  and  is  summarized  in  Table  1.  The 
known  station  and  event  geometry  is  shown  in  Figure  1,  in  addition  to  the  calculated  likelihood  function  0  -  based 
on  baekazimuth  constraints  only  -  as  a  function  of  X  and  y  with  /()  and  V  fixed  at  their  optimum  values.  Figure  2 
is  similar  to  Figure  1  exeept  that  the  arrival-time  based  likelihood  function,  O,  is  plotted  for  comparison. 


Tabic  I.  Summary  of  the  Utah  Test  and  Training  Range  explosion  ground  truth  and  observations  used  in  this 
article  to  illustrate  the  IBL  algorithm. 


Ground  truth  data 

Location 

41  13152°N,  1 12.89577°W 

Origin  Time 

2007/08/27  20:43:12.0 

Observations 

BGU 

Arrival -time 

2007/08/27  20:44:27 

Baekazimuth 

30.96° 

EPU 

Arrival-time 

2007/08/27  20:45:35 

Baekazimuth 

237.80° 

NOQ 

Arrival -time 

2007/08/27  20:47:17 

Baekazimuth 

304.22° 
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e  log(e) 


0  100  200  0  100  200 
EW(km)  EW(km) 


Figure  1.  Solutions  obtained  for  the  Utah  Test  and  Training  Range  explosion  summarized  in  Tabic  1. 

Left:  Plot  of  the  backazimuth  function  0  over  x  and  y  with  t0  and  V  Fixed  at  their  optimum 
values.  Right:  Plot  of  the  natural  logarithm  of  0,  again  with  /0  and  v  fixed  at  their  optimum  values. 


«  log(  <x) 


EW(km)  EW(km) 


Figure  2.  Left:  Plot  of  the  arrival  time  function  over  x  and  y  with  /0  and  V  fixed  at  their  correct  values. 
Right:  Plot  of  the  natural  logarithm  of  O,  again  with  tQ  and  v  fixed  at  their  optimum  values. 

The  posterior  likelihood  function  (obtained  using  Bayes’  theorem,  outlined  above)  is  shown  in  Figure  3.  It  is 
interesting  to  note  that,  whereas  in  Figures  I  and  2  (for  backazimuth  and  travel-time  constraints  individually)  we 
obtain  ellipses,  in  Figure  3  (once  these  separate  constraints  have  been  effectively  combined)  we  observe  a  circular 
region  that  is  better  constrained.  This  result  highlights  the  value  of  incorporating  both  backazimuth  and  arrival  time 
constraints  into  the  solution. 
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Figure  3.  Left:  The  posterior  likelihood  function  of  the  variables  ,V0  and  V0.  For  this  array  geometry,  the 
likelihood  contours  produced  by  the  backazimuths  run  perpendicular  to  those  produced  by  the 
arrival  times.  The  result,  seen  in  this  Figure,  is  a  well-eonstrained,  circular  source  region.  Right:  The 
natural  logarithm  of  the  same  posterior  likelihood  function. 


Collection  of  Long-Term  Ground-Truth  Events  from  a  Limestone  Mine  in  Korea 

Surface  blasting  activities  at  industrial  mines  are  good  examples  of  point  sources  that  produce  impulsive  infrasonie 
signals.  Well-defined  blasting  information  can  be  utilized  to  study  the  characteristics  of  infrasound  propagation  in 
the  atmosphere  and  to  verify  infrasonie  source  location  methods.  In  particular,  long-term  collection  of  ground-truth 
events  from  a  mine,  conducting  daily  blasting,  provides  the  basis  for  assessing  spatial  and  temporal  variations  of 
infrasound  propagation  depending  on  seasons. 

Figure  4(a)  shows  the  locations  of  several  seismic  stations  (blaek  triangles)  and  two  seismo-acoustie  arrays 
(CHNAR  and  ULDAR)  in  South  Korea.  The  inner  box  illustrates  the  distribution  of  seismic  events  near  the  cast 
coast  of  South  Korea  from  2005  to  2006.  The  Korea  Earthquake  Monitoring  System  (KEMS)  reviewed  more  than 
900  events  most  with  seismie  magnitudes  of  less  that  A/L  2.0.  Of  these  events,  those  followed  by  distinct  infrasonie 
signals  detected  by  infrasound  arrays  were  discriminated  as  surface  explosions,  not  earthquakes,  and  are  marked  by 
yellow  circles  in  the  figure.  A  field  survey  identified  a  large  open-pit  limestone  mine  (red  dot  in  the  figure)  in  the 
scattered  distribution  of  these  seismic  events.  This  mine  blasts  several  tons  of  explosives  nearly  every  working  day 
and  can  be  considered  as  an  infrasonie  bell  since  many  infrasonie  signals  from  the  mine  have  been  recorded  at  the 
seismo-acoustie  arrays  in  Korea. 

Two  temporary  seismo-acoustie  stations  were  deployed  approximately  -300  m  apart  from  the  blasting  point  during 
two  years  starting  from  April  2007  in  order  to  colleet  detailed  ground-truth  information  at  the  mine. 
Three-eomponent  seismometers  and  acoustic  sensors  were  used  to  document  the  sources.  The  acoustic  sensor  was 
connected  to  a  wind-noise  reducer.  Seismie  and  acoustic  data  were  digitized  at  100  samples  per  second. 
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Figure  4.  (a)  Locations  of  the  mine  (red  dot),  seismic  stations  (black  triangles),  and  infrasound  arrays 

(CHNAR,  ULDAR).  Inner  box  shows  distribution  of  seismic  events  around  the  mine  in  2005-2006. 
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Figure  4.  (b)  Infrasonic  signals  and  PMCC  results  for  a  delay-fired  shot  at  the  mine,  recorded  at  CHNAR, 
ULDAR,  and  near-field  seismo-acoustic  station. 


Preliminary  Analysis  of  Infrasound  Propagation  Using  the  Ground  Truth  Information 

During  the  two  years  of  sensor  operation,  the  mine  conducted  1066  blasts,  averaging  more  than  one  blast  per  a  day. 
All  blasts  were  delay-fired  shots  and  the  total  explosives  used  were  up  to  1 1.2  tons  per  blast.  The  temporary  stations 
recorded  the  near-field  seismic  and  acoustic  pressure  changes  from  these  blasts,  thereby  providing  data  to  constrain 
the  exact  origin  time,  location  and  blast  information. 
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To  quantify  the  seasonal  variation  in  infrasound  propagation  we  used  infrasound  array  observations  at  two  different 
observation  locations,  CHNAR  and  ULDAR.  Each  array  is  located  at  a  similar  distance  from  the  mine,  181  km  and 
169  km,  respectively.  However,  the  propagation  path  environments  are  quite  different.  CHNAR  is  located  to  the 
NW  from  the  mine,  producing  a  continental  path  over  land.  The  propagation  direction  for  ULDAR  is  nearly  opposite 
in  direction  to  CHNAR,  eastward  from  the  mine,  and  across  open  ocean.  Figure  4(b)  displays  recordings  of  a 
delay-fired  shot,  recorded  at  the  seismo-acoustic  arrays  (CHNAR  and  ULDAR),  and  the  temporary  station. 

Using  the  complete  dataset,  we  estimated  the  origin  time  of  each  blast  from  ncar-field  seismic  recording,  bottom 
waveform  in  the  figure.  We  applied  the  PMCC  method  to  the  distant  array  data  to  detect  the  blast-associated 
infrasonic  signals.  Wave  parameters  were  estimated  and  included  arrival  time,  azimuth,  apparent  velocity,  and 
amplitude  of  detected  signals  at  each  array.  Horizontal  propagation  velocity  of  infrasound  signals,  also  called 
celerity,  is  used  not  only  in  infrasonic  location,  but  also  infrasonic  phase  identification.  In  general,  this  velocity  can 
be  estimated  by  ray-tracing  through  a  realistic  atmospheric  velocity  model.  In  the  absence  of  a  realistic  atmospheric 
model,  the  velocity  can  be  estimated  using  ground  truth  events  with  known  origin  times  such  as  in  this  study.  The 
horizontal  velocities  for  all  the  blasts  in  the  ground  truth  dataset  were  estimated  based  on  ground-truth  source  time 
and  infrasound  arrival  times  at  both  arrays.  Figure  5  displays  the  variation  of  infrasound  travel  times  from  the  mine 
to  two  arrays  over  a  period  of  two  years. 

Based  on  this  analysis  wc  conclude: 

•  Although  propagation  distances  are  similar,  travel  times  of  the  first  arrival  (main  phase)  are  significantly  different 
along  the  paths  to  the  two  arrays.  For  CHNAR,  horizontal  velocities  range  from  0.26  to  0.29  km/s  corresponding 
to  stratospheric  returns.  On  the  other  hand,  infrasound  signals  arrive  at  ULDAR  with  fast  horizontal  velocities, 
-0.34  km/s  corresponding  to  tropospheric  phases. 

•  For  CHNAR,  travel  times  vary  systematically  with  the  seasons,  arriving  up  to  68  seconds  faster  in  summertime. 
Thus,  the  horizontal  velocity  showed  cyclic  variations  with  velocities  near  -0.26  km/s  in  the  winter-spring  time, 
and  increasing  to  -0.29  km/s  in  the  summertime,  and  decreasing  to  -0.27  km/s  in  the  autumn.  This  variation  in 
velocity  indicates  that  travel  times  of  infrasound  signals  can  be  strongly  dependent  on  seasonal  ambient  sound 
velocity  structure  in  the  atmosphere. 

•  The  travel  times  to  ULDAR  show  little  seasonal  variation.  There  is  some  diurnal  variation  of  horizontal  velocity. 
These  observations  may  be  a  result  of  the  dominant  effect  of  the  ocean  on  temperatures  at  the  ocean-atmospherc 
boundary. 

•  These  path  variations  contribute  to  signal  detectability  as  well.  Signals  at  CHNAR  have  the  highest  detectability  in 
the  summertime,  but  many  signals  were  not  detected  during  autumn-winter-spring  seasons,  even  though  strong 
source  signals  were  generated  at  those  times.  Compared  with  CHNAR,  detectibility  of  ULDAR  was  low  only  in 
summer,  but  stable  during  the  other  seasons.  Thus  the  seasonal  detectability  of  the  two  arrays  was  opposite  to  each 
other,  as  shown  in  the  figure. 
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Figure  5.  The  variation  of  travel  time  of  infrasound  arrivals  from  the  mine  to  CHNAR  and  ULDAR  during 
the  two  years  of  the  ground-truth  dataset.  Note  that  ULDAR  was  installed  in  December  2007. 
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CONCLUSIONS  AND  RECOMMENDATIONS 


We  have  developed  the  framework  for  implementing  a  statistical  infrasound  event  locator.  At  present,  the  algorithm 
is  implemented  in  a  Cartesian  coordinate  system,  suitable  for  local  to  regional  distances,  awaiting  straightforward 
conversion  to  a  spherical  coordinate  system.  The  algorithm  provides  a  thorough  characterization  of  the  source  region 
defined  by  likelihood  contours  around  the  estimated  source  location.  We  are  in  the  process  of  assigning  probability 
values  to  these  contours  by  numerically  integrating  the  posterior  likelihood  function  (which  has  a  total  integral  equal 
to  1).  In  an  example  case,  we  show  that  the  combination  of  both  arrival  time  and  backazimuth  constraints  results  in 
an  improved  solution  over  using  either  constraint  separately. 

The  preliminary  analysis  of  the  ground  truth  data  set  indicates  that  travel  times  of  infrasound  strongly  depend  on 
seasons  as  well  as  path  environments.  In  particular,  infrasonic  waves  to  CHNAR  propagated  as  guided  waves 
between  the  ground  and  stratosphere.  Infrasonic  waves  to  ULDAR  also  propagated  as  guided  waves,  but  observed 
horizontal  velocities  indicate  that  wave  ducts  were  formed  at  relatively  lower  heights  in  the  troposphere,  even 
though  propagation  distances  to  these  arrays  are  similar.  The  data  also  show  that  detectability  is  dependent  on 
season.  Detectability  in  summer  is  higher  than  spring-winter-autumn  seasons  in  the  direction  to  CHNAR  with 
ULDAR  showing  reverse  detectibility.  These  results  indicate  that  infrasound  propagation  is  related  to  path 
environment  and  seasonal  wind  in  and  around  the  Korean  Peninsula.  We  interpret  the  time  variation  in  detectability 
as  partially  controlled  by  prevailing  periodic  winds,  northwesterly  in  the  winter  and  southeasterly  in  the  summer.  To 
investigate  these  characteristics  in  detail,  further  research  is  required  to  compare  the  measured  horizontal  velocities 
with  predicted  ones  from  atmospheric  modeling.  These  results  also  illustrate  that  seasonal  and  path  dependent 
celerities  may  improve  infrasound  location  procedures. 
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